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Abstract 

For the reliable analysis and modelling of astrophysical, laser-produced and fusion plasmas, atomic 
data are required for a number of parameters, including energy levels, radiative rates and electron impact 
excitation rates. Such data are desired for a range of elements (H to W) and their many ions. However, 
measurements of atomic data, mainly for radiative and excitation rates, are not feasible for many species 
and therefore calculations are needed. For some ions (such as of C, Fe and Kr) there are a variety of 
calculations available in the literature, but often they significantly differ from one another. Therefore, 
there is a great demand from the user community to have data 'assessed' for accuracy so that they can 
be confidently applied to the modelling of plasmas. In this paper we highlight the difficulties in assessing 
atomic data and offer some solutions for improving the accuracy of calculated results. 
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1 INTRODUCTION 



Atomic data, including energy levels, radiative rates, and electron impact excitation rates, are required for the 
modelling of a variety of plasmas, such as astrophysical, laser-produced and fusion. These data are also needed 
for determining plasma parameters, such as temperature, density and chemical abundance. Since experimental 
values for most of the above atomic parameters are neither available nor can be easily measured, one has to 
depend on theory. For ions of heavy elements (mainly Z > 28) there is a scarcity of atomic data, but for lighter 
ions often multiple calculations are available in the literature. However, significant differences are frequently 
observed among different yet comparable calculations, for almost all required atomic parameters - see, for 
example pQ. For this reason there is a great demand from the user community to have data 'assessed' for 
accuracy, so that they may be confidently applied. However, there are many difficulties in assessing atomic 
data. In this paper we discuss some of these difficulties and offer possible solutions for improving the accuracy 
of calculated data. 

2 ATOMIC PARAMETERS 

The atomic parameters most commonly required for the modelling of plasmas are the following. 

2.1 Energy Levels 

The transition energy E^- and wavelength A^- are related by the following simple relationship: 

E ij = E j - Ei = hvij = hc/Xij (1) 

where i and j are the lower and upper levels of a transition, h is Planck constant, c is the velocity of light in 
vacuum, and i/y is the transition frequency. Atomic physicists generally refer to transition energies, whereas 
astronomers are more comfortable using wavelength terminology. For many ions experimental values have 
been compiled and critically assessed by the National Institute of Standards and Technology (NIST). Their 
assessed data are regularly published in refereed journals and are also freely available on their website: 
http://www.nist.gov/prnl/data/asd7cf3 However, often their data are incomplete (i.e. there are many 
missing levels, see for example, He-like ions), or the levels of a state have non-degenerate energies. In such 
circumstances there is no choice except to use theoretical energy values. 

2.2 Radiative Rates 

Radiative rates (A- values), also known as Einstein A-coefficients or transition probabilities, are related to the 
absorption oscillator strengths (f- values) by the following relationship: 

/« = 8^ A ^ A * = L49 X 10 " 16 4(^M)^ ( 2 ) 
where m and e are the electron mass and charge, respectively, and uii and ljj are the statistical weights of 
the lower (i) and upper (j) levels, respectively. The f- values are dimensionless and A- values are in units of 
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sec -1 . Generally, A- values for electric dipole (El) transitions are dominant but contributions from other type 
of transitions, namely magnetic dipole (Ml), electric quadrupole (E2) and magnetic quadrupole (M2), can be 
appreciable, and therefore are required in a complete plasma model. Measurements of A- values are sparse 
and therefore these are often determined theoretically. Unlike energy levels, the A- values are not rigorously 
assessed by the NIST and the corresponding data available on their website are very limited. 

2.3 Lifetimes 

The lifetime r for a level j is defined as follows: 
1 



Since this is a measurable parameter, it provides a check on the accuracy of the calculations for A- values. 
However, measurements are confined only to a few levels of a limited set of ions. Theoretically, lifetime 
calculations include mainly the contribution of the El transitions, but those from other types of transitions 
may be significant - see, for example, Aggarwal and Keenan 

2.4 Electron Impact Excitation Collision Strengths 

The collision strength is related to the better-known quantity collision cross section (tr) as follows: 

ttij(E) = kfaiO-ijiTval) (4) 

where k\ is the colliding energy of the electron, u}{ is the statistical weight of the lower (i) level, and ao is 
the Bohr radius. Since is a dimensionless quantity, intercomparisons among various calculations become 
straightforward. As for A- values, measurements of cr or are very limited and therefore hardly provide 
any comparisons with theory. Furthermore, calculations and/or measurements for f2 at a few energies are 
not sufficient, as the thresholds energy region is often dominated by numerous closed-channel (Feshbach) 
resonances - see Figs. 6-9 of pQ. Therefore, values of il need to be calculated in a fine energy mesh in order 
to accurately account for their contribution. Furthermore, in a hot plasma electrons have a wide distribution 
of velocities, and therefore values of f2 are generally averaged over a Maxwellian distribution to determine the 
effective collision strengths as follows: 

T(T e ) = / n(E)eM-E J /kT e )d(E J /kT e ), (5) 

where k is Boltzmann constant, T e is the electron temperature in K, and Ej is the electron energy with respect 
to the final (excited) state. Once the value of T is known the corresponding results for the excitation q(i,j) and 
de-excitation q(j,i) rates can be easily obtained - see Eqs. (9-10) of pQ. Effective collision strengths do not vary 
strongly with changing electron temperature, and therefore it is easier to fit them to a polynomial function of 
T e , as in the CHIANTI database at http://www.damtp.cam.ac.uk/user/astro/chianti/. The contribution 
of resonances may enhance the values of T over those of the background values of collision strengths (fis), 
especially for the forbidden transitions, by up to an order of magnitude (or even more) depending on the 
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transition and/or the temperature. Similarly, values of need to be calculated over a wide energy range 
(above thresholds) in order to obtain convergence of the integral in Eq. (5). 



2.5 Line Intensity Ratio 



The intensity of an emission line can be expressed as: 



AjiNjNAjNAhv. 
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where Nj is the relative population of level j, Na,z is the relative ionic abundance of ion with charge Z 
of element with atomic number A, is the relative (with respect to hydrogen) chemical abundance, Njj e 
is the relative chemical abundance of helium, n is the total number density of hydrogen and helium nuclei 
(in cm -3 ), and L is the path length through the line emitting region. However, this equation applies to 
astrophysical plasmas whereas for laboratory and fusion plasmas some parameters, such as L and N# e are not 
required. Nevertheles, calculations for the intensity of a single emission line requires many parameters, which 
subsequently add to the uncertainties. Therefore, the ratio of two lines i — > j and m — > n of an ion, i.e. 



eliminates many parameters as it then depends only on the wavelengths, A- values, and populations of the 
upper levels for which the above noted (and some other parameters, such as ionisation, recombination, and 
photoexcitation cross sections) are required. This is the main reason that a set of lines is often used as plasma 
diagnostics. If the two lines of an ion have common upper levels (i.e. j = n), then R normally depends only 
on the wavelengths and A- values, except e.g. in the presence of opacity. However, such lines are not very 
useful as diagnostics because they are independent of the density and/or temperature of the plasma. Those 
which are useful vary with either density or temperature. 



For calculations of the above atomic parameters a variety of structure and scattering codes are available. Some 
of the most commonly used atomic structure codes to generate energy levels and A- values are: Configura- 
tion Interaction Version 3 (CIV3: [4 ), General- purpose Relativistic Atomic Structure Package (GRASP: [5]), 
Superstructure (SS: 6 ), AutoStructure (AS: [7]), Multi- Configuration Hartree-Fock (MCHF: [5]), Many-Body 
Perturbation Theory (MBPT: [9 ), Hartree-Fock Relativistic (HFR: 10 J, Flexible Atomic Code (FAC: [TT]). 
and Hebrew University Lawrence Livermore Atomic Code (HULLAC: j!2j). Some of these codes have been 
published (such as CIV3, SS and GRASP), some are available on request (HULLAC), and some are on the 
web (FAC: http://sprg.ssl.berkeley.edu/~mfgu/fac/). Similarly, the most commonly used scattering 
methods are the R- matrix and distorted-wave (DW). The associated codes are the standard i?-matrix (RM: 
|13|). which also incorporates the one-body relativistic operators, and the fully relativistic version, i.e. the 
Dirac Atomic R-matrix Code (DARC: http://web.am.qub.ac.uk/DARC/). There are many versions of the 



R = 




(7) 



3 CODES FOR CALCULATIONS OF ATOMIC DATA 



5 



DW method, but the most commonly used are those of University College London [TJ], FAC and HULL AC. 
Furthermore, FAC and HULLAC are self-sufficient codes providing data for a range of atomic parameters, 
whereas the RM and DARC require input wavefunctions from CIV3 (although SS and AS can also be used) 
and GRASP, respectively. Finally, it must be stressed that some codes are comparatively easier to use than 
others (particularly FAC), but with some practice it is not too difficult to obtain the required results from the 
other codes. However, the desired level of accuracy for the results is difficult to achieve. It is even more difficult 
to determine the accuracy, as we will discuss later. Finally, almost all codes are under constant development 
and therefore the latest (and the most suitable) version may differ from those published or made available 
earlier. 

4 CHOICES FOR A MODEL CALCULATION 

For any model calculation a choice has to be made for a number of parameters some of which are listed below. 

1. Number of states/levels 

2. Configuration interaction (CI) 

3. Inclusion/exclusion of pseudo states 

4. Inclusion/exclusion of relativistic effects (RE) 

5. Energy and temperature ranges 

6. Number of partial waves 

7. Top-up 

8. Inclusion/exclusion of resonances 

9. Elimination of pseudo/spurious resonances 
10. Inclusion/exclusion of radiation damping 

The inclusion/exclusion of the above parameters is generally crucial for the accuracy of the generated 
parameters. However, which to include/exclude and at what level/s depends on the requirements, such as: 
(i) where the data are to be applied, (ii) what level of accuracy is desired or acceptable, (iii) what codes are 
available to the worker/s and/or with what codes they are familiar with, and most importantly (iv) what 
computational resources and time-frame are available. Depending on these requirements a choice has to be 
made before beginning a calculation, particularly for the determination of T. This is because some calculations 
are so computationally challenging and time consuming that they cannot be easily repeated. For the same 
reason, i.e. to economise on computational resources, some of the above noted factors have to be compromised, 
which often leaves scope for improvements in subsequent work. However all the above parameters may not 
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be important for a particular species. For example, CI is generally more important for lighter ions (around 
Z < 20) - see, for example, Aggarwal [15] , but may sometimes be equally important for heavier ions - see, for 
example, Aggarwal et al [16] for energy levels of Ni XIII to Ni XVI. In some cases, particularly for H-like ions, 
the inclusion of a large number of configurations is either not sufficient or not feasible, and therefore pseudo 
states need to be included [17] to achieve the desired accuracy in the derivation of energy levels and subsequent 
parameters. Often, extensive CI is included for the determination of energy levels and radiative rates, but 
restrictions are imposed for the further calculation of collisional data due to the limitation of computational 
resources. However, this approach leads to the generation of pseudo resonances [18] which need to be smoothed 
over a wide energy range to avoid the overestimation of T. 

Another important contribution to the determination of energy levels is from the relativistic operators, 
namely mass correction, spin-orbit interaction, spin-other-orbit interaction, spin-spin interaction, Breit inter- 
action, Darwin term, and quantum electrodynamics (QED) effects. Generally, the heavier the ion the more 
important is the contribution of these operators - see for example [19] . However, sometimes their significance is 
appreciable even for lighter ions, as shown by Aggarwal et al |20[ . Particularly important is their contribution 
in splitting the fine-structure levels of a term - see, for example, Hamada et al [5T] for energy levels of H-likc 
ions. Finally, for most of the lighter ions (Z < 20) inclusion of one-body relativistic operators, in a Breit-Pauli 
approximation, is sufficient as noted by [22], but for heavier ions fully relativistic calculations are preferable 

mm- 

For collisional calculations the choices are crucial for the inclusion of: (i) the number of states/levels, (ii) 
the number of partial waves, and (iii) the energy range up to which values of Q are calculated. This is because 
electron-ion scattering, especially at low energies, can be considered as a two-step process, i.e. a temporary 
capture of the colliding electron by the target, followed by autoionization [13] : 

e~ + A, -> (A-)* -> e~ + A 3 (8) 

where the asterix indicates the (N+l) resonance states. It is the description of such processes which is explicitly 
included in i?-matrix, but is (generally) ignored in the DW method. The contribution of resonances may 
dominate the determination of T at lower temperatures, particularly when they arise close to the thresholds 
- see, for example, Figs. 1-3 of Aggarwal [23] for transitions among the fine-structure levels of the ls 2 2s 2 2p 2 
ground configuration of Ne V. For these transitions, resonances have enhanced the values of T by over an order 
of magnitude at a temperature of 10 4 K. However, sometimes the contribution of resonances can be significant 
even at temperatures as high as 10 6 K, as demonstrated by Aggarwal and Keenan [24] for transitions in Mo 
XXXIV and by Aggarwal et al [35] for Gd XXXVII. This is particularly true when the energy difference 
between any two levels is very large. For example, the 2s 2 2p 5 2 Si/ 2 and 2s 2 2p 4 3s 4 Ps/2 levels of Mo XXXIV 
are separated by over 155 Ryd (see Table 1 of [24]). and the entire energy region is dominated by resonances, 
as shown in Fig. 1 of Aggarwal and Keenan [2"4"j. 

Certain types of radiative rates from recombination resonances increase with ionization stage, while the au- 
toionizing rates remain relatively constant [26 . As a result, in many highly ionized species radiative rates from 
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resonances can become comparable to the autoionizing rates, i.e. the dielectronic recombination competes with 
electron impact excitation, and hence can cause a significant reduction in the resonance contributions. This 
effect is called radiation damping, and can be significant for some transitions, as discussed and demonstrated 
by [26], [22] for He-like ions and by Ballance and Griffin [28] for W XXXXVII. Nevertheless, a majority of 
transitions are (generally) not affected by radiation damping, particularly at the high temperatures at which 
data are required for plasma modelling, as demonstrated by [27] j [2S] an d discussed by Aggarwal and Keenan 

mm- 



5 ASSESSMENT OF ATOMIC DATA 

The assessment of atomic data is a tedious and never ending job, as newer and newer data keep appearing 
in the literature. However, it requires expertise and experience to assess any data. Furthermore, it is not 
always easy to get the assessed data published, which inhibits the workers to undertake the job. Nevertheless, 
several institutions have undertaken the work of data assessment in the past. For example, in the 1970s, 
Los Alamos National Laboratory (USA) collected and assessed atomic data, Queen's University Belfast (UK) 
did the similar work in the 1980s, and in the 1990s National Institute for Fusion Science (Japan) was quite 
active in storing, assessing and publishing assessed data. At present, the National Institute of Standards and 
Technology (USA) is the only one active in storing, assessing and disseminating data, but this is restricted to 
energy levels and is mostly from measurements, although they do also compile some A- values. As there has 
always been a great demand from the user community to have assessed data, several databases have appeared 
in the past to fill the void, but only partially - see section 8 below. However, in the recent past with the efforts 
of the IAEA, and keeping in view the future requirement of the fusion community (particularly ITER), the 
National Fusion Research Institute (South Korea) has voluntarily offered its services for data assessment. 

One of the major difficulties in assessing atomic data, particularly for radiative rates and collisional pa- 
rameters (0 and T), is the lack of measurements, because these are not easy to perform. Furthermore, even 
if laboratory measurements are available for a few transitions and at a few energies for f2, these are not 
very helpful to assess T as data are required over a wide range of energies. Therefore, for most of the ions 
theoretical data are the only choice. Where a few calculations are already available, as for a few Fe ions, it 
is comparatively easier to assess the accuracy of the data. The major difficulty for assessment arises when 
there is a single calculation - see section 7 below for details. Some of the calculations are so computationally 
demanding and time consuming that they require months (if not years) of work, even with the best resources 
available - see for example [28] for W XXXXVII. Since it is not feasible to repeat the calculations, if an error 
occurs in the basic input data, it cannot be easily corrected. However, in a majority of cases large errors 
can be avoided simply by making extensive comparisons. If no previous calculation is available with which to 
compare then it is advisable to perform a parallel calculation with the FAC code, as it is freely available, easy 
to implement and quick to run. The alternative is to do some modelling with the help of databases and to 
compare the end results, such as line intensities or their ratios, with the observations. However, this approach 
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is not always helpful as observations may be scarce, taken at different times with separate instruments and 
cover different wavelength ranges, or lines may be blended. Nevertheless, we discuss below some of the errors 
commonly noticed in atomic data and offer the reasons for likely discrepancies between calculations. 

6 TYPES OF ERRORS 

Generally, large discrepancies between any two sets of data are due to three types of errors, namely (i) inherent 
error/s in the code, (ii) non (rigorous) assessment of data, and (iii) the (unjustified) approximations made in 
the calculations. The first type is difficult to detect but easier to correct. Often errors are noted with the 
repeated use of code/s for a variety of ions or following comparisons of results with other similar calculations, 
as noticed for He- like pQ , [55J , [30] and Li- like [31], [35] ions. This is the main reason that actively used codes 
are often under continuous development. The second type of error appears mainly when author/s either make 
unrealistic assumptions or assess the accuracy of their data based on expectation rather than rigorous tests 
and comparisons. As examples, see [33J,[33] for Ni XIX and [35] for Ni XI, which demonstrate the inadequate 
inclusion of partial waves and the energy ranges. Finally, in all calculations some approximations have to 
be made, as already discussed in section 3. This (generally) leaves scope for further improvement /s in the 
calculated data, and hence is a continuous process, but large discrepancies are rarely noted. 

Before performing calculations for an ion it is important to keep the application in mind, because different 
plasmas have different requirements. For example, many photoionized astrophysical plasmas (such as H II 
regions and planetary nebulae) have electron temperatures below 50,000 K, and hence the position of near- 
threshold resonances are very important because of their dominant contribution at low temperatures - see 
Figs. 1-3 of Aggarwal [53] and a recent work by Palay et al [3j5]. On the other hand, collisionally dominated 
plasmas (such as solar and stellar coronae) have temperatures ~10 6 K, which means that the position of 
near-threshold resonances may not be crucial, but calculations for f2 need to be performed over a large energy 
range in order to achieve convergence of the integral in Eq. (5). Similarly, a choice of a larger model of the ion 
(i.e. enlarging the number of levels included in a calculation) may significantly alter (improve) the subsequent 
calculations of T, because of the inclusion of resonances arising from the additional levels - see particularly 
the recent work }37j-|40| on Fe ions. Finally, laser-produced and fusion plasmas may require the calculations 
of T up to ~10 8 K, depending on the ion. As stated earlier, the contribution of resonances is more appreciable 
in calculations of T at comparatively lower temperatures, but their contribution at higher temperatures can 
also be significant as noted in the cases of Kr XXXII [41], Mo XXXIV [24] and Gd XXXVII [25]. Ions of some 
elements, such as C, Si and Fe, have applications in a variety of plasmas. However, a calculation performed 
with one particular application in mind (e.g. for astrophysical plasmas) may not be directly applicable to 
another type of plasma. Doing this may lead to large errors in the analysis. 
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7 DISCREPANCIES IN ATOMIC DATA 

In principle, calculations performed for any atomic parameter by using any method/code should give similar 
results (if not the same), provided the model sizes (and other associated parameters) are comparable. Unfor- 
tunately however, that is often not the case. Here we discuss the kind of discrepancies frequently found for 
different parameters and their likely causes. 

7.1 Energy Levels 

As an example, energy levels of the 2s2p 6 3£ configurations of Ni XIX from the CIV3 code are lower by 
up to 1.5 Ryd than those calculated by the GRASP and FAC codes, or experimentally compiled by NIST 



(http://www.nist.gov/pml/data/asd.cfm) - see Table 1 of Aggarwal and Keenan [32] ■ The main reason 
for such a large discrepancy is that Hibbert et al [33] focused their attention only on the lowest 27 levels of 
Ni XIX, but reported results for up to 37 levels. On the other hand, energy levels for Ni XVII from the CIV3 
calculations [44] differ by over 2 Ryd from other theoretical and experimental results, as confirmed by our 
calculations not only from the GRASP and FAC codes, but also from the CIV3 - see Table 5 of Aggarwal et al 
[45] . In this case the differences in energy levels are not understandable and there must have been some error 
in the calculations of [44] . Similarly, energies of [27] from the AS code for levels of He-like ions are higher by 
up to 2 Ryd [19] , depending on the ion, because two-body relativistic operators were not included in the code. 
Thus there can be several reason/s for the discrepancy in the reported energy levels. 

Sometimes there is a problem in the assignment of the orderings of energy levels, as seen in Table 1 of Gupta 
et al [46] for Al-like ions. The 2 D^ 2 5/2 energy levels of the 3p 3 and 3s3p( 3 P°)3d configurations are interchanged 
in the MBPT calculations of Safronova et al [47] for a series of ions as confirmed by the experimental results as 
well as our calculations from three independent codes, namely CIV3, FAC and GRASP. The most likely reason 
for the interchange is that Safronova et al performed calculations for a large number of ions but investigated 
the orderings for only a few (lighter) ions. In many instances the orderings of levels are not the same for all 
ions in a series, and hence the discrepancy. Generally, listings at the NIST website are helpful in assigning the 
level orderings, but occasionally there are differences with the theory as noted in the cases of Kr XXXI and 
Kr XXXII [48]. 

In all calculations, the designation of an energy level is (mainly) determined on the basis of the strength 
of its corresponding eigenvector. However, often and particularly for those ions whose levels extensively mix 
with different configurations, such as Fe XVI [49] . it is not always possible to assign a unique (unambiguous) 
designation for a level, because a single eigenvector may dominate for several levels. In such cases, the best one 
can state is that a particular level has such a J value of parity even or odd, but the corresponding configuration 
from which it comes remains ambiguous, and subject to the interpretation of the individual author/s and/or 
user/s. 
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7.2 Radiative Rates 

If the same level of CI (and relativistic operators) are included in a calculation, the f- or A- values should 
agree within ^20%, irrespective of the method or the computer code employed, and this is true particularly 
for strong (f > 0.01) transitions. For weaker transitions, the discrepancies among different calculations can 
be large, often by an order of magnitude or even more depending on the transition. This is because weaker 
transitions are more sensitive to different levels of CI and/or their energies (Ai?,j), and for the same reason 
their length and velocity forms also often differ quite significantly. On the other hand, the strong transitions 
are generally more stable in magnitude, and their length and velocity forms also agree closely, which gives 
an indication of the accuracy of a calculation. However, discrepancies for strong transitions can also be large 
among different calculations as seen in Table 5 of Aggarwal and Keenan [42] for Ni XIX. The differences in f- 
values obtained from the GRASP and FAC codes are up to 50% with those from CIV3 [43], particularly for 
transitions involving levels 28 and higher, due to the corresponding differences in energy levels as discussed 
above. Similarly, f- values from the CIV3 code [50] differ from the GRASP and FAC results, by up to an 
order of magnitude, for many transitions of Fe IX as may be seen in Table 5 of Aggarwal et al [51] . These 
differences arise in spite of the fact that the energy levels of [50] are in close agreement with the experimental 
values, because of the 'adjustment' of the Hamiltonian (known as 'fine-tuning'). It is difficult to fully explain 
such large differences in A- values for so many transitions, but one can speculate on the reason/s. In most 
calculations from the CIV3 code, after a preliminary survey all those levels/configurations whose eigenvectors 
(mixing coefficients) are very small (say < 0.01) are removed from the final calculations in order to economise 
on computational effort, and this process affects the weaker transitions more than the stronger ones, because 
of the additive or destructive effect of the components. The choice of an appropriate cut-off level below which 
the levels/configurations are removed depends on the authors/s and/or on the size of the calculations, and we 
believe this is the main reason for the large differences in f- values discussed above for transitions in Fe IX. 
Therefore, from this (and many other similar) examples we may conclude that the process of fine-tuning may 
make theoretical energy levels more accurate in magnitude, but not necessarily the subsequent calculations of 
radiative and collision rates. Similarly, the orderings of the levels may remain uncorrected in the absence of 
experimental energies for a larger number of levels, if not all. 

Differences in A- values for several (particularly weaker) transitions of Li-like ions are up to three orders of 
magnitude, as demonstrated by Aggarwal and Keenan [31] . These types of large discrepancies are difficult to 
detect without performing independent calculations using different code/s. While the application of inaccurate 
data may affect the modelling of plasmas, it may not always be possible to perform calculations with different 
codes. However, caution may be exercised by performing several tests with the same code. 

7.3 Collision Strengths and Effective Collision Strengths 

Adopting the same level of complexity, different scattering methods and/or codes should provide comparable 
values of f2, for a majority of transitions. However, differences among several comparable calculations are often 
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abnormally striking as we discuss here with some examples. As shown in Fig. 4 of Aggarwal and Keenan |52) . 
values of from the DW calculations [S3] are significantly lower than those from darc for several forbidden 
transitions of Fe XVI, and the discrepancy between the two calculations increases with increasing energy. This 
is because Cornille et al |53) included only a limited number of partial waves in their calculations with the 
assumption this would be sufficient for the convergence of f2, which is not the case as may be seen from Figs. 
1-3 of Aggarwal and Keenan [S3] . Eissner et al [SI] corrected this limitation in their calculations from the R- 
matrix code, but their subsequent results of T are still underestimated for several transitions over a wide range 
of temperature as seen in Fig. 12 of Aggarwal and Keenan [S3]. This is because their calculations involved 
only the lowest 21 levels of Fe XVI, while the larger calculation of Aggarwal and Keenan [S3] with 39 levels 
also included resonances arising from the higher levels. However, even the results of Aggarwal and Keenan 
may be further improved in a similar way by the inclusion of yet higher levels of the n > 6 configurations of 
Fe XVI. Improvements over a previous calculation is a continuous process, particularly with the increasing 
availability of computational resources. 

Fe XVI is a moderately heavy ion for which the contribution of rclativistic effects is important, but not 
dominant, in the calculation of atomic parameters. However, the semi-relativistic i?-matrix calculations of 
Bautista |55) for T differ from the fully relativistic results from DARC by up to an order of magnitude, for 
several transitions, as noted in Table 2 of Aggarwal and Keenan [56^. This is because of the sudden shift in the 
background values of fl (£Ib) by the inclusion of relativistic effects, as shown in Fig. 3 of Bautista [55]. The 
shifts in Qb are upwards, downwards, and random, and are not realistic, but happened because of an error in 
the adopted version of the code which the author was not aware of at the time of the calculations. Most of 
the atomic codes currently in use are (generally) always under continuous development, as already stated and 
is illustrated by the several versions of the FAC code available on the web. An error in a code is often difficult 
to spot and harder to speculate in advance, and can only be noted and corrected after its prolonged use for a 
variety of calculations, and mainly by comparisons with other similar but independent calculations. 

Generally, differences of up to a factor of two or so in values of 51 between any two calculations, at a 
few energies, do not considerably affect the subsequent determination of T. Nevertheless, discrepancies in T 
are sometimes up to three orders of magnitude, as recently demonstrated by Aggarwal and Keenan [T] for 
transitions of several He-like and Li-like ions [31], [32]. Particularly affected are those transitions which belong 
to the degenerate levels of a state (often referred to as 'elastic' transitions because of their very small energy 
differences). Furthermore, these differences in T are not confined to the lower or higher values of temperature, 
but persist over the entire range as can be seen in Figs. 11 and 12 of Aggarwal and Keenan [ST]. Such large 
discrepancies in T reported by Whiteford et al [27] happened, yet again, due an error in the adopted version 
of the i?-matrix code, and has now been corrected [58j . However, without the availability of independent 
calculations with the DARC code, it was not possible to know the errors in the reported data for several He-like 
[27] and Li-like [59] ions. 

Calculations of (and subsequently T) for allowed transitions among degenerate levels of a state, such as 
2si/2 - 2p!/2,3/2 in H-like ions, are difficult to perform, because of their very slow convergence |21] with partial 
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waves. Additionally, a slight variation in their AE^ may easily lead to an under/overestimation in f2 and 
subsequently of T, as discussed and demonstrated in detail by Hamada et al |21) . 

Even without an error in the adopted version of a code, values of T may significantly disagree among 
calculations by different authors as may be noted from Figs. 12-18 of Aggarwal and Keenan [60] for transitions 
of O IV. Some of these discrepancies arise due to the non-convergence of either f2 (because of the limited number 
of partial waves included) or the integral in Eq. (5), because of the limited range of energy included, as clearly 
demonstrated in Figs. 2, 3 and 7 of Aggarwal and Keenan [35] for transitions in Ni XI. Resonances close to the 
thresholds can easily move left or right of a threshold due to differences in energy levels, and their positions 
sometimes considerably affect the calculations of T, particularly at very low temperatures (~ 10 3 -10 4 K). 
This is the main reason for the differences in the behaviour of the variation of T with T e as noted above for 
transitions of O IV. 

8 SOURCES OF DATA 

Atomic data in the literature are often spread over a wide range of journals and this makes the task of even 
identifying data difficult for a user, apart from the problems of the assessment of the results as discussed 
above in sections 6 and 7. Additionally, due to a significant increase in computing and storage power during 
the past decade or so, a typical calculation generates so much data that no journal can publish these in 
their entirely, although some do provide a significant amount in their electronic versions. To overcome this 
difficulty, a few websites store a significant amount of data, for a variety of parameters and for a large range 
of ions. The most common and widely used websites, especially for collisional data, are: (i) CHIANTI: 
http://www.damtp.cam.ac.uk/user/astro/chianti/, (ii) ADAS: http://open.adas.ac.uk/, (iii) APAP: 
http://amdpp.phys.strath.ac.uk/ and (iv) CFADC: http://www-cfadc.phy.ornl.gov/. Apart from the 
numerical atomic data, these websites also provide a wide range of computer programs for the applications of 
the data, and hence are extensively employed by those who generate, assess, and/or apply the atomic data 
to the modelling of plasmas. Therefore, particularly for a user who is not an expert in atomic data, these 
websites and repositories of data are very helpful. However, websites also have some disadvantages. Data 
may be incomplete, are not assessed by independent experts, and are often unpublished. The other problem 
is that revising/updating the website is a continuous process, and since this is staff intensive, there can be 
a considerable delay before the (more accurate) published data are incorporated into the data repository. 
Furthermore, in rare cases data pertaining to the same calculation may be different on different websites, as 
recently highlighted by Aggarwal and Keenan [I] for transitions of Fe XXV. As a consequence, the data on a 
website may not always be the best available in the literature. Therefore, for an active researcher who wants 
to use the latest and/or the best available atomic data, there is no choice except to search for the data himself. 
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9 ADVICE 



In this paper we have discussed a few examples (although there are many more) of discrepancies for different 
atomic parameters. There can be several reasons for such discrepancies among different calculations. However, 
we would like to state that a calculation can only be considered to be the 'best' available until a better one can 
be performed, as there is always scope for improvement. For example, by (i) considering a larger model of an 
ion, (ii) including more CI in the generation of wavefunctions, (iii) widening the range of partial waves and/or 
energy, (iv) resolving resonances in a narrower energy mesh, (v) including radiation damping, (vi) including 
two-body relativistic operators (if not already done) , and (vii) including the effects of higher lying ionization 
channels through pseudostates, which are particularly important for lighter H-like ions, such as He II and Li 
III. In conclusion, improving a calculation is always a continuous process and depends on the availability of the 
computer (and staff) resources, but large discrepancies (errors), particularly for radiative rates and collisional 
data, can be avoided (or considerably reduced) by a few simple precautions as stated below for the benefit of 
the producers, assessors, and users of atomic data. 

9.1 Producers 

1. Make as many comparisons as possible between different calculations for a variety of transitions, such 
as: allowed, forbidden, semi-forbidden, weak, and strong. 

2. In the case of large discrepancies, try to understand and explain these without making assumptions. 

3. Report results for collision strengths (Q), at least for a few transitions and at a few energies, so that a 
relationship with the subsequent results for T can be formed. Most of the discrepancies in values of T 
can only be understood by knowing the corresponding differences in Q. 

9.2 Assessors 

1. This is a difficult task to perform (although none is trivial), especially when it is not easy to assess even 
one's own work. 

2. Assess what methods and assumptions have been used in a calculation. 

3. Follow some basic guidelines, such as: behaviour of a transition, adequacy of the J/L and E ranges, 
inclusion (exclusion) of resonances, relativistic effects, etc. 

9.3 Users 

1. The best situation is when only one set of data is available. However, this is often not the case. 

2. If two or more data sets are available, and the authors do not fully and convincingly justify the improve- 
ments made, then use both (or more) sets of data and make your own assessment. However, remember 
that the latest available calculation may not always be the best. 
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3. In case of doubt and/or suspicion, contact the authors. 
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